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Abstract — The performance of second order statis- 
tics (SOS) based semi-blind channel estimation in 
long-code DS-CDMA systems is analyzed. The co- 
variance matrix of SOS estimates is obtained in the 
large system limit, and is used to analyze the large- 
sample performance of two SOS based semi-blind 
channel estimation algorithms. A notion of blind es- 
timation efficiency is also defined and is examined via 
simulation results. 

I. Introduction 

In practical wireless communication systems, the transmit- 
ted signals usually experience fading, which either attenuates 
the received power or causes dispersion. Usually, the chan- 
nel state information is unknown to both the transmitter and 
the receiver in many practical applications, thus necessitating 
channel estimation at the frontend of the coherent receiver. 
In this paper, we focus on the long-code direct-sequence code 
division multiple-access (DS-CDMA) systems and frequency- 
selective fading channels. 

A variety of channel estimation algorithms have been de- 
rived, based primarily on one of two aspects of random ob- 
servations, namely, the distribution and the moments. When 
the distribution of the received signal conditioned on the chan- 
nel state information is known, the maximum likelihood (ML) 
criterion can be applied to yield asymptotically optimal per- 
formance. Such ML channel estimation algorithms are suit- 
able for training-symbol-based systems, in which a fraction of 
the transmitted symbols is known to both the transmitter and 
receiver. However, the channel state information hidden in the 
information symbols is ignored, thus decreasing the spectral 
efficiency. 

When only the information symbols are available (usually 
called blind channel estimation), the moment based channel 
estimation algorithms require substantially less computational 
cost than other methods, thus being feasible for practical ap- 
plications. Thus, over the past two decades, a large number 
of blind channel estimation algorithms [T][7]|ll | |13 ) have been 
developed using moment estimation, particularly the second 
order statistics (SOS). Typically, based on SOS estimation, 
the subspace method [Q and moment matching can be 
applied. The subspace technique is suitable only for stationary 
channels, e.g. short-code CDMA systems, since it requires the 
signal subspace to be time-invariant. However, in many prac- 
tical CDMA systems, long codes are employed, thus making 
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the overall channel non-stationary. Therefore, only the mo- 
ment matching technique can be applied directly to long-code 
CDMA systems. When training symbols are available, both 
the training and information symbols can be used jointly to 
yield better estimates (usually called semi-blind channel esti- 
mation) . 

In this paper, we adopt the moment-matching-based algo- 
rithm of 11 II for estimating the SOS of channel coefficients of 
each user. The following two SOS based semi-blind channel 
estimation algorithms are considered. 

• Moment-matching-based estimation. The asymp- 
totically optimal moment based estimator [Sj is applied 
to the first order moment, estimated from the train- 
ing symbols, and the SOS, estimated from the informa- 
tion symbols, thus resulting in moment-matching-based 
semi-blind channel estimation under some assumptions. 

• Subspace-based estimation. The subspace spanned 
by the channel coefficient vector is obtained from the 
eigenstructure of SOS estimation. The channel esti- 
mates are obtained by projecting the training symbols 
onto the subspace followed by linear combining. It 
should be noted that the subspace here is different from 
that mentioned above, which is spanned by the signal. 

The remainder of this paper is organized as follows. A 
discrete model of the frequency selective fading channels is 
explained in Section II. The performance of SOS estimation 
and semi-blind channel estimation is analyzed in Section III 
and Section IV, respectively. Simulation results and conclu- 
sions are given in Section V and Section VI, respectively. The 
details can be found in the journal version of this paper 

II. System Model 

We consider a synchronous uplink DS-CDMA system operat- 
ing over frequency selective fading channels of order P (i.e, 
P is the delay spread in chip intervals). Let K denote the 
number of active users, N the spreading gain and (3 = ^ the 
system load. We model the frequency selective fading chan- 
nels with discrete finite-impulse-response (FIR) filters. The 
z-transform of the channel response of user k is thus given by 

p-i 

M*) = £s*(pV> (1) 

where {gk(p)} P =o,...,p-i are the corresponding independent 
and identically distributed (i.i.d.) (with respect to both k 
and p) channel coefficients with variance -p. We adopt the 
block fading model, in which the channel coefficients, namely 
{gk(p)}, remain unchanged during a coherent block, and we 



denote the coherence time by M, measured in channel symbol 
periods. 

With the assumption P <C N, we can ignore the effects of 
intersymbol interference (ISI). Thus, the chip matched filter 
output at the l-th chip period in the m-th symbol period can 
be written as 



r(mN + l)=^2 x k {m)hi m \l) + n{mN + I), 



(2) 



where {xk(m)} denotes the complex channel symbols, which 
are mutually independent with respect to k and m. We as- 
sume that Xk(m) is a complex random variable satisfying 
E{\x k (m)\ 2 } = 1 and E {x|(m)} = 0, which hold, for ex- 
ample, for quadratic phase shift keying (QPSK) modulation. 
{n(mN + 1)} is additive white complex Gaussian noise, which 
satisfies E{\n(mN + l)\ 2 } = a 2 x . {h^ (I)} denote the con- 
volution of the spreading code and the channel coefficients: 



(l)*9k(l), 



(3) 



where s^ (I) is the Z-th chip of the original random spreading 
codes of user k in symbol period m, which takes values -H= 
and — equiprobably and independently with respect to k, 
m and I. 

When the channel coefficients are unknown to the receiver, 
channel estimation is necessary for coherent reception. For 
simplicity, we assume that M% training symbols, which are 
known to both transmitter and receiver, are transmitted at 
the beginning of each coherent block, and the proportion of 
training symbols is denoted by a 



M t 
M 



III. Performance of SOS Estimation 

When discussing SOS estimation in this section, for notational 
simplicity, we assume that all M channel symbols during the 
coherent block are unknown to the receiver. When analyz- 
ing the performance of semi-blind channel estimation, we can 
simply replace M with M — Mt in the conclusions for SOS 
estimation. It should be noted that our analysis on the per- 
formance of SOS estimation is based on the large system limit, 
namely K, N — > oo while keeping /3 and P constant. 

A SOS Estimation Algorithm 

We adopt the algorithm proposed in jllj to estimate the 
channel coefficients in long-code CDMA systems. On defin- 
ing 2 r(m) = (r(mN + P), r((m + 1)N)) T as the signal un- 
affected by ISI during the symbol period m, we can rewrite 
in vector form, which is given by 



r(m) = ^2 C£ m) gfcXfcfm) + n(m) 



(4) 



where g k 



(g k (0),..., 9k (P~l)y 

(m) 



tifc(m) 



(n(mN + P),...,n((m + l)N)y and C k m > is the truncated 
Sylvester matrix, which is given by 



,(m) 



>(P) 



(P- 



(1) 

»(2) 



>(N - P + l) J 



^^Note that cr^ is the noise variance, normalized to represent the 
inverse signal-to-noise ratio. 

2 Superscript T denotes transposition and superscript H denotes 
conjugate transposition. 



We define a cost function for moment matching, which is 
given by 



J : 



J2 (E{H r (m)}-H r (m)) 



(•») 



where R r (m) = r(m)r(m) H and || ■ \\f denotes the Frobenius 
norm. The optimal SOS is obtained by minimizing the cost 
function J, 

d = arg min J, 

d 

where d = (df , d k )" , d k = vec (gkSk ) and vec(X) de- 
notes a vector formed by stacking the columns of matrix X 
into one column vector. An explicit expression for d is ob- 
tained in which is given by 

d = T- x y, 

3 

M 



(6) 



where 



- P 2 KxP 2 K 



-L J2 Q T (m)Q(r 



A I 



1 

yp^Kxi = J7 Yl Q T ( m ) Yec (r(m)r H (m))-CT 2 Q T (m)vec(I) 



Q(m) ( 



JV-P+l) 2 xP 2 K 



(Q 1 (m),Q 2 (m), Q K (m)) , 



and 



Qk(m)(N-p+i) 2 xp 2 — C k 



(m) 



c 



(m) 



where ® denotes the Kronecker product. In practical imple- 
mentations, d can be obtained by solving @ iteratively. 

B Ident inability and Consistency 

On analyzing the elements of T, we obtain the following 
lemma. 

Lemma III.l T converges elementwise to I almost surely as 
M — » oo and N — > oo. 

In it is shown that the channel is identifiable and that 
SOS based channel estimation is consistent if T is nonsingu- 
lar. We cannot draw the conclusion that T is nonsingular as 
AI — > oo and N — > oo merely from Lemma 1111.11 since the 
corresponding convergence is elementwise and not in norm. 
However, we can show the following proposition, which states 
a sufficient condition for the non-singularity of i?{T}. 



Proposition III. 2 Suppose K and N are such that KP < 
N - P + 1. Then E{T} is nonsingular and ||£{T} _1 - i)^ 
is finite 4 . 

As K and N become sufficiently large, the condition that 
KP < N - P + 1 is equivalent to /3 < |j. Since T con- 
verges to -E{T} almost surely as M — ■> oo, T is nonsingu- 
lar with large probability when M is sufficiently large. Since 
||i5{T} _1 — l|| is finite as K, N — + oo, most of the elements 
of _E{T} -1 converge to the corresponding elements of I. Thus 
it is reasonable to approximate T _1 with I when K and N 
are sufficiently large. For the case f3 > -p, numerical simu- 
lation indicates that the non-singularity of T still typically 
holds. Thus, in the following analysis of SOS estimation er- 
ror, we assume that T = I, which will be validated later by 
simulation results. 



3 The dimensions of the matrices are labeled as subscripts. 
*||^4nxn||l = maxi<j<„ 



C SOS Estimation Error 

The SOS estimation error Sdk = dfc — dfc is due to the multiple- 
access interference (MAI) and noise. The following proposi- 
tion gives an explicit expression for the covariance matrix of 
Sdk, where the assumptions can be shown to hold in elemen- 
twise convergence as K, N — > oo. 

Proposition III. 3 Assume 

K 

^ Qy("^)dydyQ i j(7n) = /3 2 l (N _p +1 - ) 2 x(N _ P+l - ) 2. 



IV. Performance of Semi-blind Channel 
Estimation 

In this section, we consider two approaches to semi-blind 
channel estimation, namely moment-matching-based and 
subspace-based algorithms. In semi-blind channel estimation, 
Mt training symbols are available (thus, the number of infor- 
mation symbols is M — Mt) and the corresponding received 
training signal, which is represented by a iVMt-vector, is given 
by 



K 



(8) 



(Cl m) ® n(m)) sld (C ( k m) ® n(m)) = (3o 2 n 



fc = 2 



where, similarly to Cj(fc) is the Sylvester matrix con- 
structed by 



£ (n*(m) ® C< m) ) g fe gf (n*(m) ® Ci m) ) = /3a 



i fc = ( (s«x fe (l) 



Then we have, for k — 1, K , as K,N — > oo, 

2/3 



By multiplying j^C? (k) on both sides of ©, the MAI 
for user can be completely cancelled, as Mt — > oo since 
-j^Cf (k)Ct(k) — > I. Therefore, a sufficient statistic for gfc 

is given by gfc = -j^pl (fc) r t- Thus, the equivalent training 
signal for user k is given by 



where 



gfc = gfc + n fc , 



(9) 



|(gfe)|-^j| + |(Sfc)mod(i,p)| • tf* = J 
(gfc)r i i (gfc)^j i otherwise, if mod(i, P) — mod(j,P) 

(g fc )mod(i,p) (g fc )mo<i(j,f>) ' 
0, otherwise 



otherwise, if — ^-j^j 



An interesting observation on (|7J is that the SOS estima- 
tion errors are mutually correlated and the covariance matrix 
is dependent on the realization of channel coefficients. How- 
ever, when (3 is moderate or large, the correlation is weak, 
compared with the error variance. We can apply the asymp- 
totic results in J7J to finite systems as approximations. 

For the average SOS estimation variance of all users, we 
have 



1 K 



KP 2 



from which we can see that the performance of SOS estimation 
is determined by the system load and noise variance. 

It should be noted that it is not sufficient to obtain 
E {<5dfc<5dj? } for the performance analysis of semi-blind chan- 
nel estimation, since the channel coefficients are complex. 
E {<5dfc<5djT} is also necessary for the analysis in the next sec- 
tion. However, it is easy to check that, for all 1 < i,j < P 2 , 



£ { 5d ^}), 3 = ^( EL ), 



(mod(j,P)-i)p+[f] 



where rife = -^Cf (k)nt is additive white complex Gaussian 

CT 2 

noise with variance jg-. As Mt — > oo, (fife) ; become mutually 
independent with respect to k and I. 

Since both {gfc} and {dfc} across different users are weakly 
correlated when M is sufficiently large , we can carry out the 
channel estimation for each user separately, thereby reducing 
considerably the computational complexity at cost of marginal 
performance loss. For notational simplicity, we drop the sub- 
scripts of user index throughout this section. The elements in 
g and d are denoted by g p , p = 1, P, and d p , p = 1, P 2 , 
respectively. 

A Moment Matching 

A.l Moment-matching-based estimator 

For applying the asymptotically optimal estimator in 8 
and the corresponding analysis, it is necessary to discuss 
the moment based channel estimation in R. For any com- 
plex vector v, we denote the corresponding real vector by 

v = (jft (v) T , SJ (v) T ^j , where 5R (v) and Si (v) represent the 
real and imaginary parts of v. 

It should be noted that there are only F ^ +1 ' > free vari- 
ables in d and P 2 free variables in d (note that there are P 
real elements in d, corresponding to the diagonal elements in 
gg T ) since d is the vectorization of a Hermite matrix. We 
denote the real vector of these free variables by d / , where the 
subscript / means free. 



Then, it is further easy to check that (E {5d k 8dl}) . . = 0, 
i = 1,...,P 2 , except for \-p\ = mod(i, P) which means that 
(dfc)^ are the diagonal elements of gg H , and thus are real. 
Therefore, we can conclude that the real and imaginary parts 
of the SOS estimation error are uncorrelated and have identi- 
cal variances. 



On defining the observation vector z 



(g J , d J ) T and Sz = z— z, the asymptotically optimal estimator 
using estimates of the first and second moments |S] is obtained 
by minimizing a cost function, which is given by 



Jopt (g) = (z (g) 



(2(g) 



(10) 



where T,$ z = E {<5z(5z T } and the function z maps g to z. 
It is easy to check that 



Ex, 



i: 



I2 







1 



M-M t ' 



where the covariance matrix E^dy can be obtained from Prop. 

ITTO 

However, due to the conclusion of Prop. MI. 31 T,sa f is de- 
pendent on the realization of g, thus being unknown to the 
estimator. Therefore, the optimal estimator is infeasible for 
practical applications. Since the cross correlation of different 
elements in d/ is small for moderate or large /3, we can assume 



that E« 



a d I, thus resulting in the cost function given by 



(11) 



J (g) = w vec I gg I — d + (1 — w 



where the weighting factor w — 



(l-q)g-^ 



Since minimiz- 



ing 1111 is equivalent to obtaining g, which optimally matches 
the first and second moment estimates in z, we call it moment 
matching based channel estimation. 

A. 2 Performance analysis 

Practically, the cost function 1111 can be minimized with it- 
erative optimization methods. For theoretical analysis, the 
optimal g can also be obtained by taking derivatives of the 
cost function mil with respect to g, resulting in 2P equations 
denoted by Fi = 0, i — 1, 2P. These equations determine 
a mapping ^Sf from the observation z to g in a neighborhood 
of z. Thus, by applying the implicit function theorem we 
can obtain the Jacobian matrix of vE" at z, which is given by 



dz 



dg 



dz 



(12) 



, F2p) , provided that |p- is non-singular. 



where F = (fi, *—^ T - A " A +w 81 

Therefore, for the channel estimation error <5g = g — g, we 
have the following proposition, whose proof is essentially the 
same as that of Theorem 3.16 in 8 and is omitted in this 
paper. 

Proposition IV. 1 In moment-matching-based semi-blind 
channel estimation, if \fMSz converges weakly to a random 
vector with zero mean and a covariance matrix Y*sz as M — > 
00, then MSg converges weakly to a random vector with zero 
mean and covariance matrix given by 



E,, a - I — Sj,- ' 



dz 



V dz 



(13) 



that the Jacobian matrix ^S. 

az 



A. 3 Performance bound 



is nonsingular. 



Using Lemma 3.1 in ||, the asymptotic performance of SOS 
based blind channel estimation can be lower bounded by that 
of the asymptotically optimal estimator, which is given by 



where 



/ I2 



OA, 



\ 



(15) 



B Subspace-based Approach 
B.l Subspace-based estimator 

Another methodology for semi-blind channel estimation is to 
make use of the subspace estimated from the SOS estimates. 
In this section, we adopt a simple subspace algorithm, which 
can be carried out in the following three steps: 

1. The one-dimensional subspace spanned by g is obtained 
from the unit eigenvector u of matrix vec -1 (d) 5 cor- 



responding to its largest eigenvalue. 

The channel estimate from the training symbols, namely 
g, is projected onto the subspace, thus decreasing the 
noise power and obtaining a tentative channel estimate, 
which is given by 



H - 

u gu. 



(16) 



3. The final channel estimate is obtained by combining g 
and g linearly: 

g = wg + (1 - w)g, 

where ui is a weighting factor, which will be optimized 
later. 

This algorithm can also be regarded as tackling the phase 
ambiguity incurred by SOS based blind channel estimation by 
making use of the training symbols. 

B.2 Perturbation on signal subspace 

On denoting the acute angle between the subspace spanned 

( Is^ul 

by g and the unit eigenvector u by 9 — arccos [ M „ M | , we 



can obtain an explicit asymptotic expression for E^sm 2 ^, 
which is given in the following Proposition. 

Proposition IV.2 As K,N,M ~* 00, 

o'l = Af£;|sin 2 e} 



(P - 1) ((2/3 + t|g|| 2 ) *l + K) 2 + P 2 + 2 -f) 



(17) 



C Performance of channel estimation 



On defining Su = 



S "iig|| g " ' **1^ can ^ e rewritten as 



|u g| 



+ Su 



; + n 



+ Su 



where ii is defined in @. It is easy to check that 5u H Su a 
sin 2 8, which becomes tight as 8 — * 0. We assume that the 
different elements in Su are mutually independent, and the 
real and imaginary parts of each element are also mutually 



independent with identical variances. These assumptions can 
be validated by simulation. 

By ignoring the higher order terms and applying |u H g| — » 
||g|| 2 , we can obtain 

o ME\s g H Sg\ 



■ 2 "sll 2 (i + £)<rg , (i-^) 2 (P-i) 



az (7 



(1 - a)P Pa 

where the optimal u) is given by 



+ 7^, (18) 
Pol 



l|g£(l+j>| (P-l)o 
1 — cm a. 



(19) 



When uj = 0, a g = — , which equals the performance of 
training-symbol-based estimation. Thus, the subspace-based 
channel estimation with optimal ui surely attains better per- 
formance than training symbol based estimation. When the 

2 o- 2 

subspace estimation is perfect and lu = 1, a g — — p, which 
means that noise level is lowered by a factor of -p. Therefore, 
the last term in 1181 is the performance with perfect SOS 
estimation and the first term in 1181 represent the penalty 
incurred by imperfect SOS estimation. 

V. Simulation results 

A Moment-matching-based Channel Es- 
timation 

We define the blind estimation efficiency as 



(20) 



Intuitively, r\ quantifies how many training symbols each infor- 
mation symbol is equivalent to for the purpose of semi-blind 
channel estimation. Figure 1 shows the efficiency r\ obtained 
from Prop. II V. II versus /3 and a\ when P — 3 and a — 0.2. 
We observe that r\ decreases monotonically with (3 since larger 
implies more MAI for SOS estimation. An interesting ob- 
servation is that the efficiency is small for both high and low 
noise levels and achieves a maximum (around 0.3) for mod- 
erate noise power (around a\ = 1). An intuitive explanation 
is that, when a\ is large, the SOS estimation performance 
is poor since the SOS error variance contains a term propor- 
tional to (ofj) ; when a\ is small, the SOS estimation error is 
dominated by the MAI, thus achieving much lower efficiency 
than training-symbol-based estimation. 

Figure 2 shows the efficiency versus P and a when a\ =0.5 
and /3 = 0.5. An important observation from this figure is 
that, for large a and P, the efficiency r\ < 0, which means that 
it is harmful to incorporate the information symbols into the 
channel estimation. This arises from the suboptimality of the 
channel estimator in 11 U since the exact covariance matrix 
of the SOS estimation error is unknown and the estimated 
weighting factor w is also suboptimal. Thus, the unreliable 
SOS estimation may have negative impact on the reliable es- 
timation from the training symbols, when a is large. When 
a is small, the efficiency achieves its maximum for small P; 
however, the efficiency does not change much with respect to 
P. 



From the simulation results of Fig. 1 and Fig. 2, we can 
draw a conclusion that moment-matching-based channel esti- 
mation is suitable for systems with small system load, small 
numbers of training symbols, small channel order and moder- 
ate noise power. 

B Subspace-based Channel Estimation 

The blind estimation efficiency is shown in Fig. 3 (versus f3 
and o~\ when a — 0.1 and P = 3) and in Fig. 4 (versus 
a and P when f3 = 0.5 and o~\ = 0.5). Similar to that of 
moment-matching-based channel estimation, the optimal effi- 
ciency is attained at moderate noise levels (around = 1), 
and the efficiency decreases with respect to /3. However, unlike 
moment-matching-based channel estimation, the efficiency of 
subspace-based channel estimation increases with a and P. 
Thus, subspace-based channel estimation is suitable for sys- 
tems with large channel order and large numbers of train- 
ing symbols, which are opposite to the desirable conditions 
for moment-matching-based channel estimation. Also, we see 
that the blind estimation efficiency of subspace-based channel 
estimation is always positive, indicating that this technique 
is more robust than moment-matching-based channel estima- 
tion. 

VI. Conclusions 

In this paper, we have analyzed the performance of SOS based 
semi-blind channel estimation in long-code CDMA systems. 
The main results include the following: 

• An explicit expression for the covariance matrix of the 
SOS estimation error has been obtained in the large 
system limit within some assumptions. 

• Expressions for the performance of two types of semi- 
blind channel estimation have been obtained. Particu- 
larly, we have obtained an asymptotic expression char- 
acterizing the perturbation of the eigenvectors of covari- 
ance matrices, which can be applied to other problems 
as well. 

• The blind estimation efficiency has been obtained from 
simulation results, which show that the SOS based semi- 
blind channel estimation attains high efficiency for sys- 
tems operating in the moderate noise region. Moment- 
matching-based estimation is suitable for systems with 
small channel order and small numbers of training sym- 
bols, while subspace-based estimation achieves good 
performance in systems with large channel order and 
large numbers of training symbols. 
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Figure 3: Blind channel estimation efficiency versus noise 
power and system load in subspace-based channel estima- 
tion 
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Figure 1 : Blind channel estimation efficiency versus noise 
power and system load in moment-matching-based chan- 
nel estimation 




Figure 2: Blind channel estimation efficiency versus chan- 
nel order and training symbol proportion in moment- 
matching-based channel estimation 



